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Abstract 


The excellent mechanical properties of laminated composite structures make them 
prime candidates for wide variety of applications in aerospace, mechanical and other 
branches of engineering. The enormous design flexibility of advanced composites is ob- 
tained at the cost of large number of design parameters. Due to complexity of the systems 
and lack of complete design based informations, designers tent to be conservative in their 
design. Furthermore, any new design is extensively evaluated experimentally until it 
achieves the necessary reliability, performance and safety. However, the experimental 
evaluation of composite structures are costly and time consuming. Consequently, it is 
extremely useful if a full-scale structure can be replaced by a similar scaled-down model 
which is much easier to work with. Furthermore, a dramatic reduction in cost and time 
can be achieved, if available experimental data of a specific structure can be used to 
predict the behavior of a group of similar systems. 

This study investigates problems associated with the design of scaled models. Such 
study is important since it provides the necessary scaling laws, and the factors which 
affect the accuracy of the scale models. 

Similitude theory is employed to develop the necessary similarity conditions(scaling 
laws). Scaling laws provide relationship between a full-scale structure and its scale model, 
and can be used to extrapolate the experimental data of a small, inexpensive, and testable 
model into design information for a large prototype. Due to large number of design param- 
eters, the identification of the principal scaling laws by conventional method( dimensional 
analysis) is tedious. Similitude theory based on governing equations of the structural 
system is more direct and simpler in execution. The difficulty of making completely sim- 
ilar scale models often leads to accept certain type of distortion from exact duplication 
of the prototype (partial similarity). Both complete and partial similarity are discussed. 
The procedure consists of systematically observing the effect of each parameter and corre- 
sponding scaling laws. Then acceptable intervals and limitations for these parameters and 
scaling laws are discussed. In each case, a set of valid scaling factors and corresponding 
response scaling laws that accurately predict the response of prototypes from experimen- 
tal models is introduced. The examples used include rectangular laminated plates under 
destabilizing loads, applied individually, vibrational characteristics of same plates, as well 
as cylindrical bending of beam-plates. 
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DISCUSSION 


The importance of employing small scale models in designing advanced composite 
structures has been gaining momentum in recent years. With a view to better under- 
standing the applicability of these models in designing laminated composite structures, 
an analytical investigation was undertaken to assess the feasibility of their use. Em- 
ployment of similitude theory to establish similarity among structural systems can save 
considerable expense and time, provided the proper scaling laws are found and validated. 

Before small scale models can be used, the technical barriers that must be overcome 
are: 


• What are the proper scale factors. 

• What is the effect of these scale factors(scale effect). 

In this study the limitation and acceptable interval of all parameters and corresponding 
scale factors are investigated. In the present studies, the material behavior was assumed 
to be linearly elastic. Therefore, scale effects are not present. 


Completed Tasks^ 

An analytical investigation has been conducted in order to establish the applicability 
of similitude theory to laminated rectangular plates. Particular emphasis is placed on the 
case of free vibration and buckling of plates under uniaxial compressive and shear loads. 
Angle ply, cross ply and quasi-isotropic configurations were chosen for investigation. 


Qurrerd Tasks 

i) Nonlinear Kinematics , Linear Constitutive Relations 

The large deflection analysis is performed on composite beam- columns subjected to 
an eccentric axial compressive load. The loads are static. The objective of this investi- 
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gation is to develop and validate the scaling laws for nonlinear kinematic behavior(large 
deformation) of simple generic structural elements. 

ii) Failure Analysis 

The laminate failure analysis is investigated for first ply failure of composite beam- 
columns. The beam-columns are subjected to static eccentric compressive loads. The 
objective of this study is to develop and validate the necessary scaling laws corresponding 
to stress analysis of laminated structures. In particular we like to predict the stress profile 
in prototype by projecting the corresponding stresses of its scale model. 

iii) Curved Configuration 

The development of scaling laws which pertain to the elastic stability response of 
laminated shells is currently being investigated. Particular emphasis is placed on the case 
of buckling of orthotropic laminated cylindrical shells under axial compressive load. 

Future Tasks 

The following research and development must be done before small scale model can 
be used in design and analysis of the composite laminated structures. 

• Complete the investigation of curved configurations. 

• Study the effect of boundary conditions. 

• Study the applicability of scaled down models to geometrically stiffened structures. 

• Study the effect of geometric imperfections. 

• Design prototype and scaled down models for experimental validation( at this time it 
is anticipated that tests will be performed at the “NASA Langley Research Center” 
labs). 
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• Structural Similitude For Laminated Structures. 

• Structural Similitude and Scaling Laws for Cross-Ply Laminated Plates. 
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STRUCTURAL SIMILITUDE FOR 
LAMINATED STRUCTURES 


G. J. Simitses and J. Rezaeepazhand 
University of Cincinnati, Cincinnati, OH 45221, U.S.A. 

Abstract— Due to special characteristics of advanced reinforced composite materials, they require 
extensive experimental evaluation. Thus, it is extremely useful to use available experimental data 
of specific structural systems to predict the behavior of alt similar systems. This study describes the 
establishment of similarity conditions between two structural systems. Similarity conditions 
provide the relationship between a scale model and us prototype, and can be used to prediet t 
behavior of the prototype by extrapolating the experimental data of the corresponding small sc 
model. Since satisfying all the similarity conditions simultaneously is in most cases unpractical, 
distorted models with partial similarity (with at least one similarity condition relaxed) are 
employed. Establishing similarity conditions, based on direct use of gaming equatio s, 
discussed and the possibility of designing distorted models is investigated. The meth °^' s de ^° , 
strated through analysis of the cylindrical bending of othotropic lammated beamplates subjected 
to transverse loads and buckling of symmetric laminated cross-ply rectangular plates subjected to 
uniaxial compression. 
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NOMENCLATURE 

a plate length 

Ajj laminate extensional stiffnesses 

b plate width 

By laminate coupling stiffnesses 

D plate flexural stiffness 

D u laminate flexural stiffnesses 
Eij Young’s moduli of elasticity 
F stiffness ratio 

h total laminate thickness 

k xx bending curvature in the laminate 

m, n number of half waves in x and y 
M cross-ply ratio 

moment resultant 
N number of layers 

stress resultant 

N x inplane load per unit width 

P total transverse load 

q transverse load intensity 

Qijf Qu lamina stiffness elements 
t ply thickness 

u y Vy w reference (midplane) surface displacements 
Xy y, z reference axes 

midplane extensional strain 
0 fiber orientation angle 

A transformation matrix 

A, scale factors 

v t Vjj Poisson’s ratios 

normal stress in the kih lamina 
m model 

p prototype 

pr. predicted 

th. theoretical 


INTRODUCTION 

Aircraft and spacecraft comprise the class of aerospace stuctures that require efficiency 
and wisdom in design, sophistication and accuracy in analysis and numerous and careful 
experimental evaluations of components and prototype, in order to achieve the necessary 
system reliability, performance and safety. 
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Preliminary and/or concept design entails the assemblage of system mission require- 
ments, system expected performance and identification of components and their 
connections as well as of manufacturing and system assembly techniques. This is 
accomplished through experience based on previous similar designs, and through t e 
possible use of models to simulate the entire system characteristics. 

Detail design is heavily dependent on information and concepts derived from the 
previous step. This information identifies critical design areas which need sophisticated 
analyses, and design and redesign procedures to achieve the expected component 
performance. This step may require several independent analysis models, which, in many 

instances, require component testing. > 

The last step in the design process, before going to production, is the verification of 
the design. This step necessitates the production of large components and prototypes in 
order to test component and system analytical predictions and verify strength and 
performance requirements under the worst loading conditions that the system is expected 
to encounter in service. 

Clearly then, full-scale testing is in many cases necessary and always very expensive. 
In the aircraft industry, in addition to full-scale tests, certification and safety necessitate 
large component static and dynamic testing. The C-141A ultimate static tests include eight 
wing tests, 17 fuselage tests and seven empennage tests (McDougal, 1987). Such tests 
are extremely difficult, time consuming and definitely absolutely necessary. Clearly, one 
should not expect that prototype testing will be totally eliminated in the aircraft industry. 

It is hoped, though, that we can reduce full-scale testing to a minimum. 

Moreover, crashworthiness aircraft testing requires full-scale tests and several drop 
tests of large components. The variables and uncertainties in crash behavior are so many 
that the information extracted from each test, although extremely valuable, is nevertheless 
small by comparison to the expense. Moreover, each test provides enough new and 
unexpected phenomena, to require new tests, specially designed to explain the new 
observations. 

In the building construction industry, when the skeleton frames are erected at the 
site, a specified sequence of erection events must be followed in order to avoid collapse. 
This was discovered through (expensive) experience, but it is not widely known. A small- 
scale testing of similar structures would definitely have been safer and less costly. 

Finally, full-scale large component testing is necessary in other industries as well. 
Ship building, automobile and railway car construction all rely heavily on testing. 

Regardless of the application, a scaled-down (by a large factor) model (scale model ) 
which closely represents the structural behavior of the full-scale system ( prototype ) can 
prove to be an extremely beneficial tool. This possible development must be based on the 
existence of certain structural parameters that control the behavior of a structural system 
when acted upon by static and/or dynamic loads. If such structural parameters exist, a 
scaled-down replica can be built, which will duplicate the response of the full-scale 
system. The two systems are then said to be structurally similar. The term, then, that best 

describes this similarity is structural similitude . 

Similarity of systems requires that the relevant system parameters be identical and 
these systems be governed by a unique set of characteristic equations. Thus, if a relation 
or equation of variables is written for a system, it is valid for all systems which are similar 
to it (Kline, 1965). Each variable in a model is proportional to the corresponding variable 
of the prototype. This ratio, which plays an essential role in predicting the relationship 
between the model and its prototype, is called the scale factor . In establishing similarity 
conditions between the mode! and prototype two procedures can be used, dimensional 

analysis and direct use of governing equations. p 

Models, as a design aid, have been used for many years, but the use of scientific 
models which are based on dimensional analysis was first discussed in a paper by Rayleigh 
(1915). Similarity conditions based on dimensional analysis have been used since Rayleigh’s 
time (Macagno, 1971), but the applicability of the theory of similitude to structural 
systems was first discussed by Goodier and Thomson (1944) and later by Goodier (1950). 
They presented a systematic procedure for establishing similarity conditions based on 
dimensional analysis. 
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In the 1950s and 1960s many interesting works were published in this area. Most of 
these authors discussed similitude theory based on dimensional analysis. Kline (1965) 
gives a perspective of the method based on both dimensional analysis and the direct 
use of the governing equations. Szucs (1980) is particularly thorough on the topic of 
similitude theory. He explains the method with emphasis on the direct use of the 
governing equations of the system. 

Due to special characteristics of advanced reinforced composite materials, they have 
been used extensively in weight efficient aerospace structures. Since reinforced composite 
components require extensive experimental evaluation, there is a growing interest in small 
scale model testing. Morton (1988) discusses the application of scaling laws for impact- 
loaded carbon-fiber composite beams. His work is based on dimensional analysis. Qian 
et al. (1990) conducted experimental studies of impact loaded composite plates, where 
the similarity conditions were obtained by considering the governing equations of the 
system. These works and many other experimental investigations have been conducted to 
characterize the size effect in material behavior for inelastic analysis. 

In recent years, due to large dimensions and unique structural design of the proposed 
space station, small scale model testing and similitude analysis have been considered as the 
only option in order to gain experimental data. Shih et al. (1987), Letchworth et al. 
(1988), Hsu et al. (1989), and McGowan et al. (1990) discussed the possibility of scale 
model testing of space station geometries especially for vibration analysis. Most of these 
studies have used complete similarity between model and prototype. 

The present study presents the applicability of small scale models, especially distorted 
models, in analysing the elastic behavior of large and complex structural systems. By 
applying similitude theory, we try to find a set of conditions between two similar 
structural systems (scaling laws). Later, these conditions can be used to design a model, 
the experimental data of which can be projected in order to predict the behavior of the 
prototype. The objectives of the investigation described herein are: 

• create necessary similarity conditions in order to design an accurate distorted 
model; 

• evaluate the derived similarity conditions analytically. 

Similarity conditions provide the relationship between model and its prototype, and 
can be used to extrapolate the experimental data of a small and less expensive model 
in order to predict the behavior of the prototype. In all of our work in this area we will 
restrict ourselves to linearly elastic material behavior. Furthermore, it is assumed that the 
laminates are free of damage (delaminations, matrix cracking, fiber breaks, etc.). 

THEORY OF SIMILITUDE 

Similitude theory is concerned with establishing necessary and sufficient conditions 
of similarity between two phenomena. Establishing similarity between systems helps to 
predict the behavior of a system from the results of investigating other systems which have 
already been investigated or can be investigated more easily than the original system. 
Similitude among systems means similarity in behavior in some specific aspects. In other 
words, knowing how a given system responds to a specific input, the response of all 
similar systems to similar input can be predicted. 

The behavior of a physical system depends on many parameters, i.e. geometry, 
material behavior, dynamic response and energy characteristic of the system. The nature 
of any system can be modeled mathematically in terms of its variables and parameters. A 
prototype and its scale model are two different systems with similar but not necessarily 
identical parameters. The necessary and sufficient conditions of similitude between 
prototype and its scale model require that the mathematical model of the scale model can 
be transformed to that of the prototype by a bi-unique mapping or vice versa (Szucs, 
1980). It means, if vectors X p and X m are the characteristic vectors of the prototype and 
model, then we can find a transformation matrix A such that: 

X„ = AX m or X m = A _1 X p . 


( 1 ) 
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The elements of vector X are all the parameters and variables of the system. A 
diagonal form of the transformation matrix A is the simplest form of transformation. 
The diagonal elements of the matrix are the scale factors of the pertinent elements of the 
characteristic vector X. 

0 0 

0 k xl 0 

o o ■ ■ ■ X xn 

where X x . - x ip /x im denotes the scale factor of x<. In general the transformation matrix is 
not diagonal. 

Since similitude theory gives many alternative ways for investigating a system, it has 
been used in areas which primarily involved many experimental investigations, such as 
fluid mechanics, aerodynamics, hydraulics and modal analysis. 

In establishing similarity conditions between the model and prototype two procedures 
can be used, dimensional analysis and direct use of governing equations. The similarity 
conditions can be established either directly from the field equations of the system or, if 
it is a new phenomenon and the mathematical model of the system is not available, 
through dimensional analysis. In the second case, all of the variables and parameters 
which affect the behavior of the system must be known. By using dimensional analysis, 
an incomplete form of the characteristic equation of the system can be formulated. This 
equation is in terms of nondimensional products of variables and parameters of the 
system. Then, similarity conditions can be established on the basis of this equation. 

In this study, we consider only direct use of the governing equations procedure. 
This method is more convenient than dimensional analysis, since the resulting similarity 
conditions are more specific. When governing equations of the system are used for estab- 
lishing similarity conditions, the relationships among variables are forced by the govern- 
ing equations of the system. 

The field equations of a system with proper boundary and initial conditions 
characterize the behavior of the system in terms of its variables and parameters. If the 
field equations of the scale model and its prototype are invariant under transformation A 
and A" 1 , then the two systems are completely similar [eqn (1)]. This transformation 
defines the scaling laws (similarity conditions) among all parameters, structural geometry 
and cause and response of the two systems. 

In order to demonstrate the applicability of the method, we consider the following 
example. Suppose we want to design a reasonable (able to test) model for a large 
rectangular plate. The plate is simply supported at all edges and loaded with a uniform 
transverse load of intensity q. Assuming uniform cross section and isotropic material, the 
governing differential equations and boundary conditions are well known (Timoshenko 
and Woinowsky-Krieger, 1959): 

d 4 w ~ d 4 w d 4 w _ 
dx 4 + dx 2 d y 1 + d y 4 D * 

and B.C. at x = 0, a 


w = 0 



and at y = 0, b 


w = 0 





( 5 ) 
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For model and prototype we may write: 


d 4 w d J w OT 

Qm 

(6) 

d4, " dx'm dyi d y x m 

" An’ 

d4w, P , , d ' w P . d * w P 

~ dx; d \y\ d y p 

_ VP 

(7) 

Dp' 


where subscripts m and p refer to model and prototype, respectively. 

By defining scale factors X it the variables of the prototype can be written as 
x ip = X x .x im . The response similarity conditions between model and prototype (complete 
similarity) are determined by substitution of the X x ,x im into the differential equation of 
the prototype and by requiring that the result be the differential equation of the model 
[eqn (5)]. 


w„ 


+ 2 


dV 


d w m / A,A q m 


±c m dyl 

Equations (5) and (8) are the same if the terms in parentheses of eqn (8) are all equal. 


Aw A w A w A^ 

A* " A^A* ~ X* y ~ X D ' 


(9) 


Now to find the scaling laws from eqn (9), we have three choices. Dividing eqn (9) by 
first term, yields: 

A, = A„ A w = ^f (10) 


Dividing eqn (9) by the second term, yields: 

X = 1 ; =22- 

A X A y , A w 

and finally dividing eqn (9) by third term: 


A A 2 A 2 

Aq A X A y 


X X 4 


(ID 


( 12 ) 


Note that all three, eqns (10)-(12), are equivalent. This means that as long as 
X x = X y , the behavioral condition that relates the response factor, X Wi to the cause scale 
factor, X qy is the same for all three cases (complete similarity). 

By applying similitude theory to a specific system, the result will be a set of 
conditions among pertinent parameters (scale factors of the parameters) of this system 
and its similar models. If all similarity conditions are satisfied, the two systems are 
completely similar. Suppose the system has m variables and similitude analysis of the 
governing equations of the system defines n relationships among m unknowns (scale 
factors of these variables). If the two systems are completely similar m - n scale factors 
can be chosen freely and the values of the other scale factors are found by using the n 
similarity conditions. The arbitrary scale factors are usually chosen based on the 
experimental facility, available material, and measurement techniques. By having the 
parameters of the prototype and scale factors, the model parameters can be calculated 
easily. Often complete similarity is difficult to achieve or even undesirable. This problem 
is usually caused by limitations on conducting the experiment. When at least one of the 
similarity conditions cannot be satisfied, partial similarity is achieved. In this case, the 
model which has some relaxation in similarity conditions is called a distorted model. 
Distorted models are more practical, since relaxation of each similarity condition 
eliminates some restrictions on the model design. These relaxations in the relationship 
between two systems cause model behavior to be different from that of the prototype. 
Understanding of these relaxations (and their effect on model behavior) can be used to 
modify the model test data so as to predict the behavior of the prototype. Since each 
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- - 

variable has different influence on the response of the system, the resulting similarity 
conditions have different influence. By understanding the effect of variables and 
similarity conditions over desired intervals, the similarity conditions which have the least 
influence can be neglected without introducing significant error (Kline, 1965). 


APPLICATIONS 



In this section, as an initial effort, similarity conditions are developed in order to 
design reasonable, distorted, scale models for orthotropic laminated beamplates and 
plates. 


a. Cylindrical bending of laminated beamplates 


v! 

a.l. Deflections. We desire to find the maximum deflection of beamplates. Beam- 
plates are subjected to transverse line loads. By assuming that the displacement functions 
are independent of y, or u = u(x), i> = 0, w = w(x) (cylindrical bending), from Ashton 
and Whitney (1970) the governing differential equations and boundary conditions are 
reduced to: 

d*w qA u (n) 

dx 4 A\\Dn — 4 i 


d 3 u B n d 4 w 
dP ~ A u dx 4 ’ 

(14) 


and the B.C.s at x = 0, a are: 



w = 0, 

(15) 

j 

, du _ d 2 w 
N xx = A n -- B n-£I= 0 . 

(16) 

— 

, _ du _ d 2 w 

( 17) 


Equation (13) can be written as: 


■r r-/ - ; ;V 

d* w 

{A u D n — 5 n) = ^11* 

(18) 


By applying similitude theory, the resulting similarity conditions, eqn (18), 

are: 

- ' •" 


(19) 

k :jj : 

or 

^u^r) u = 

(20) 



(21) 

i 

Similarly from eqns (14), (16) and (17) we have: 


3 

^-A n ^u^x — 

(22) 



(23) 


The condition depicted by eqn (23) does not give any new information, since it can 
be obtained by combining eqns (20) and (22). So, eqns (20)-(22) denote the necessary 
behavioral conditions for complete similarity between the scale model and its prototype. 

— 
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a. 1 . 1 . Parenthesis: for better understanding the restrictions of eqn (20), consider the 
definition of A mn , B mn and D mn . 


^mn S (Qmn)j 0*/ Zj— 1)» 

y* 1 
1 - v 

= - I (Q mn )j(Zj - zf-i), 

1 j- l 

A™ = ^ E (&»„);(*/ - £i). 

where Zj is the coordinate of the upper surface of the yth lamina (measured from the plate 
reference surface). Let Zj - Cjh where -0.5 < Cj < 0.5 and h is the total thickness 
0* = 0, 1, N). ( Q mn )j , the transformed stiffnesses for theyth lamina are given in terms 
of the engineering orthotropic constants and the fiber orientation angle 0. 


Q - /(#> En i £22* v i 2 > g 12 ). 


This allows us to express A mn , B mn and D mn in terms of h and functions of all Q, N 
and the stacking sequence. 

A mn = hf a (Q mn ,N ), 

B mn = A 2 A(Q™,A0. 

D mn = h l f d {Q mn ,N), 


or as scale factors: 


where: 


^ A mn ~ ^h^aiQmn > ^ 0 > 

An = XlF d (Q mn ,N), 


(24) 

(25) 

(26) 


F,= 


/(g,N) p 


i = a, 6, 


Substitution of eqns (24)-(26) into eqn (20) yields: 

A(<2n , W(&i ■ ^) = A 2 (Qn. A). (27) 

Equation (27) states that the first similarity condition, eqn (20), is independent of 
total thickness of the plate, and it is only a function of material properties, number of 
plies and stacking sequence of the model and its prototype. This condition, eqn (27), is 
satisfied if the model and prototype are made of the same material with identical N and 
the same stacking sequence of the lamina. 

Now, the accuracy of the derived behavioral similarity conditions, eqns (20)-(22) is 
evaluated analytically, in order to determine the level of confidence that can be expected 
in interpreting the data from the distorted model experiments (partial similarity). 

Consider a cross-ply laminated E-Glass/Epoxy plate composed of 96 orthotropic 
layers (0/90/0/... ) w as the prototype. We desire to find the maximum deflection of the 
prototype by extrapolating the pertinent values of a small scale model. The model has the 
same stacking sequence as the prototype but with a smaller number of layers (distorted 
model). The prototype and its scale model have the following characteristics: 

prototype (0(90/0/. . .) 96 : a = 90 in. b = 100 in. h = 0.858 in. N = 96, 

model (0/90/0/... ) 16 : a = 5.0in. b = 6.139 in. /t = 0.143 in. A=16, 

X a = 18 X b = 16.29 A* = 6 X N = 6. 


scale factors 
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Fig. 1. Theoretical and predicted maximum deflections of prototype (0/90/0/...) % when model 
(0/90/0/... ) l6 is used (a £ i| * X En = A, u * 1, A ff ■ 18, X b = X. q = 16.92, X h = X s - 6). 

For simplification we assume that model and prototype have the same material 
properties (X En = A £z2 = A„ l2 = 1) and X q = X b . By employing only the similarity con- 
dition of eqn (21) (note that X P = X x X q \ therefore the condition becomes X w X Du = X x X P ) 
the theoretical maximum deflections of the model are projected in order to predict the 
maximum deflections of the prototype. Figure 1 presents the theoretical and predicted 
maximum deflections of the prototype and corresponding theoretical values of the scale 
model. The derived scaling laws can be used with a high level of accuracy in predicting the 
prototype behavior. Note that the model was designed by employing the free scaling 
factors (partial similarity). 


a. 2. Stress analysis. For the kt h lamina, the normal stress in terms of the strains and 
curvatures (cylindrical bending) is: 

= £<?(«& + zk xx ), 

(28) 

where e® x is the extensional strain on the reference surface ( z 
change in curvature of the reference surface. 

By substituting the expressions for e ^ and 

= 0) and k xx represents the 

o$* = Qnfr. x + M - zw tXX ). 

(29) 

Applying similitude theory for the normal stress, a**. 


, . ( X u X w . A w \ 

+ u ' l ’7j- 

(30) 

The resulting similarity conditions are: 


A a (*> = Xq\vX u X x r 

(31) 

A„<*> = Aq«*>AJ,AJ 2 , 

(32) 

A a w = XqW X z X w X x , 

(33) 


where A w = A*ApA^ and X u = X w X; l X Bu X2 l u [see eqns (21) and (22)1. 

For complete similarity, eqns (3 1)— (33) give the same result. However, for the 
distorted model each similarity condition gives different results. To find which one of 
eqns (31)-(33) gives the best prediction for the prototype behavior (partial similarity), the 
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Fig. 2. Predicted and theoretical normal stress a xx distributions in various layers of the prototype 
(0 6 /90 6 /0 6 ) when (0/90/0/...) 18 is used as model. 

theoretical stresses of the model are projected with each condition and compared to 
theoretical stresses of the prototype. 

Figures 2-4 present the predicted and theoretical distributions of the normal stress c xx 
in various layers of the prototype for a cylindrical bending test. It is observed that the 
predicted stresses by eqn (33) agree very well with the theoretical results. Equation (32) 
cannot predict the behavior of the prototype accurately. Equation (31) is not a suitable 
similarity condition, since its predicted data does not match the theoretical results. The 
Figures do not include the predicted stresses using eqns (31) and (32). This is purposely 
done in order to simplify the figures. Since the results came from a three-point bending 
case, it is expected that the similarity condition based on bending normal stresses, eqn 
(33), would yield accurate stress predictions. 



Fig. 3. Predicted and theoretical normal stress g„ distributions in various layers of the prototype 
(0 6 /90 6 /0 6 ) when (0 4 /90 4 /0 4 /90 4 /0 4 ) is used as model. 
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Fig. 4. Predicted and theoretical normal stress distributions in various layers of the prototype 
(0/90/0 /... ) Ai when (0/90/0/... ) 16 is used as model. 


It is necessary to note that the stress distribution in various layers of the prototype is 
completely different from the stress distribution in the model. But, the derived similarity 
condition, eqn (33), can be used successfully to predict the stress distribution in the 
prototype. 


b. Buckling of symmetric laminated cross-ply rectangular plates 


Consider the case of cross-ply symmetric laminated plates (By — 0, D 16 — D Z6 — 
A l6 = A 2 6 = 0). The plates are subjected to inplane uniaxial compression in the x 
direction (N x ). 

The buckling loads are described only by one differential equation: 


•^ll W , < xrn + 2^12 


Wxxyy + A 


22 w jyyy 


N x w° = 0, 


(34) 


where D n = A 12 + 2D 56 . 

For simply supported plates, the boundary conditions are: 


at x = 0, a 


at y = 0, b 
The solution: 


w = 0, w° xx = 0, 

w = 0, w 0 yy = 0. 


. fmnx\ . fnny\ 
w = A mn sml ■ — 1 sml — J * 

satisfies the boundary conditions and governing differential equation if: 

R, = *’[A.(f) + 2».i(j) + D “(l) (j) 

By applying similitude theory to eqn (35): 


]2 ] 2 
A m * A n 


A 4 A 2 

A n A a 


1 _ 1 m - ] - - J - - 


(35) 


(36) 
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which yields to the following scaling laws: 


1 _ 1 m 
*N Z - *D n ^2 > 


= 


A- U 


l ^l* 


0.25 


*R 2 


1 /: \ 0.25 

_ / A P n \ 

Af| 


(37) 

(38) 

(39) 


Boundary conditions do not yield any similarity condition. Equation (37) can be 
written as: 


Ajf| ' wi ’ 


(40) 


where: 


v _ ^ 


Ajc, — 11 


2 1 2 X 2 

a D,i A 6 A m 


An. A' 


7 T £^22 * 

Conditions (38) and (39) are considered as design scaling laws and condition (40) is 
considered as response scaling law . 

b.l. Complete similarity . The necessary condition for complete similarity between 
the model and its prototype is that all the design scaling laws predict the behavior of the 
prototype with the same accuracy. In other words X Rl = A* 2 . This equality is satisfied if 
[see eqns (38) and (39)]: 

(41) 






From Tsai (1964), the bending stiffnesses D mn can be expressed as a function of the 
total number of layers, Af, the cross-ply ratio, M, and stiffness ratio, F, 

D n = [(F- l)y+ 


F>2z = [(1 ~ F)V + ^ 12 ’ 

Di 2 = 12 Qn ’ 

where F = E 12 /E n = G 22 /Q 11 an d 

1 MOV - 3)[M(N - 1) + 2QV + 1)1 


V/ = 


(IV 2 - 1)(1 + M) 3 


(1 + M) 3 

where M is the cross-ply ratio and N is the total number of plies. For the common special 
case of symmetric cross-ply laminates (0/90/0/. . .)„ . in which the laminae are all the same 
thickness and have the same material properties: 

N+ 1 


M = 


N — 1 


Substituting into eqn (41) 


r <f, - i)w,_+ ri 

■ (1 - F p )y/ p + F p 

+ 

£ 

*— < 
i 

(1 - F m )y/ m + F m 


j 2 - l 2 
a Cii “ A Ou- 


(42) 


In general, by choosing the model material and using eqn (42), the number of plies 
of model ( N m ) can be determined. Since jV m must be an integer, it is difficult to satisfy 
eqn (42), therefore partial similarity with a distorted model is pursued. 
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For simplification we assume that model and prototype have the same material 
properties (A £[( = X En = A„ p = 1). If the model and prototype have the same material 
properties then X Qu - X Qxl - 1 and F p = F m = F and eqn (42) can be simplified as: 

l(F - l)y/ p + 1][(1 - F)W P + F] = [(F - l)vf OT + 111(1 - F)^« + 

^ * (W P - 1) = Wm ‘ (V'm “ 1). 


/(/V;V m ) 


Vp ‘ ~ 1 ) 

V'm * (V'm - 0 


The numerical values of /CV p , N m ) are plotted over a large range of iV m for several 
Np in Fig. 5. From Fig. 5, it is verified that complete similarity is achieved if X N = 1. It 
is also important to notice that, as the number of plies of prototype increases, X s * 1 also 
satisfies the condition f(N pt N m ) = 1 and complete similarity is achieved. Since the 
number of plies, N p and N m are integer numbers this condition exists for large N p . 

So far we proved in the special case when X N ~ 1, condition f{N p ,N m ) = 1 is 
satisfied and complete similarity is achieved. However, there are some constraints, such 
as the geometry of the model, the model material, the number of plies and the stacking 
sequence of laminates. Since this still appears to be restrictive, we proceed with the 
determination of distorted models, for which some of these restrictions can be relaxed. 



Fig. 5. Sensitivity of complete similarity condition ^ d h ^ d v * ^£> tl f° r different 

* * complete similarity). 

b.2. Partial similarity . When at least one of the design scaling laws cannot be 
satisfied, partial similarity is achieved. In this case, since each parameter has different 
influence on the response of the system, the resulting design scaling laws have different 
influence on the accuracy of the predicted response. By understanding the effect of the 
various parameters and accuracy of the design scaling laws over desired intervals, the 
design scaling laws which have the least accurate prediction can be chosen as the “right” 
type of distortion. 

The choice of the right type of distortion is investigated as follows. In each case, all 
of the model parameters except one, are chosen to be identical to its prototype. Then, the 
effect of this relaxation for a wide range of this parameter is investigated. 
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b.2.1. Number of plies: since we assume all the laminae to have equal thickness, the 
distortion in thickness is the same as the distortion in number of plies (/t toia i = N ‘ laminae)- 
In other words, thickness is only a function of the number of plies. Consider that the 
model and prototype have the same material properties (Kevlar/Epoxy 49) with a 
different number of layers (N m ^ N p ). Figure 6 presents the per cent of discrepancy 
between theoretical and predicted values of the normalized critical load (K { ) for simply 
supported rectangular plates. In these cases, the prototype is a laminated cross-ply 
(0/90/0/... )i 01 square plate with R p = 1. The accuracy of the designed distorted models 
with the same stacking sequence but with a different number of layers [(0/90/0/...),] is 
investigated. The R m is determined by using both design scaling laws [eqns (38) and (39)]. 
Figure 6 shows that, as the number of plies of the model increases, the accuracy of the 
model increases very quickly. Models with N — 1, or 3 do not have acceptable accuracy. 
Equation (39) is the best design scaling law, especially for a model with 5 < N m < 40. For 
N m > 40 all conditions yield the same accuracy. 

The accuracy of a model with N = 13 in predicting the buckling behavior of proto- 
types with N > 101 is also investigated. It is shown that, the model with N m = 13 can 
predict the critical load of any prototype with N p = 101, ...» 500 with the same accuracy 
as prototype with N p = 101. In other words, the accuracy is independent of N p when 
N p > 101. 

This study indicates that a distorted model with a smaller number of layers can 
predict the critical load of its prototype with good accuracy. 



Fig. 6. <7o discrepancy of normalized buckling loads (/T,) when N M £ N p (distortion in N ). 

b.2.2. Material: now we consider the distortion in model material. For this purpose 
two different groups are considered: isotropic materials (which include metals and 
plastics) and fiber-reinforced composites. In all cases the prototype is considered to be an 
orthotropic laminated plate. 

The model and its prototype have the same stacking sequence, number of plies 
[(0/90/0/... ) 13 ] and aspect ratio. Figure 7 presents the per cent of discrepancy when the 
model and prototype have different material properties. For the Kevlar /Epoxy prototype, 
a Boron/Epoxy, a Boron/Polymide and most of the Graphite/Epoxys can be be used as 
the model material and vice versa. However S-Glass/Epoxy is not a good choice for 
predicting a Kevlar/Epoxy prototype and vice versa. The design scaling laws of eqn (38) 
yield the best accuracy. 
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Model Materials 

X/E 49 

B/P Mod I WRD 9371 
G/E T300 SP3I3 
G/E Mod I ERLA 4617 
B/E B(4) 5SCS 
G/E HMS 3002M 
G/E Mod I ERLA 4289 
G/E A5 35C1 
G/E Mod !! 52C6 
S-Glass/E 59G1 
G/P Mod I WRD 9371 
G/E T300 5233 

0 5 10 15 20 

% Discrepancy ol XI 

Fig. 7. <7o discrepancy of normalized buckling loads (K x ) for different composite materials. Model 
and prototype have the same stacking sequences (0/90/0/...) t3 . 



Model materials 


p: Kevlar/Epoxy 49 


Rl 


m R2 



0 5 10 15 20 25 30 35 

% Discrepancy of K1 


Fig. 8. <7o discrepancy of normalized buckling load for different isotropic materials. Model and 
prototype have the same stacking sequences (0/90/0/... )tj. 


Since plastics are used extensively for the experimental study of the behavior of the 
structures, the possibility of a plastic model or in general a model with isotropic materials 
is considered. Figure 8 presents the per cent of discrepancies for models with isotropic 
materials. The prototype is Kevlar/Epoxy plate [(0/90/0/... ), 3 ]. Almost all of the 
plastics, copper and aluminum which are used, give the same accuracy for the design 
scaling law of eqn (39) (30% discrepancy). But the design scaling law of eqn (38) yields 
better accuracy (less than 5.3%). 


DISCUSSION 

An extensive study based on analytical investigations has been conducted in order to 
establish the applicability of similitude theory to simple structural elements. Theory of 
similitude is used to design scale models for orthotropic laminated beamplates and to 
predict the behavior of the prototype, with reasonable accuracy. Similarly data of scale 
models are projected to predict prototype behavior. Even for models with different 
numbers of plies and stacking sequence of layers (distorted models) the predicted data are 
well-matched with full scale prototype data. 

By establishing similarity conditions the model parameters are specified. First, there 
is a need to verify the derived similarity conditions. The verification should be based on 
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exact analytical solution of the system. Furthermore, a validation procedure based on 
experimental data of two or three small scale models of the actual system (prototype) 
should be conducted. For example, at least two different small scale models of the proto- 
type should be designed and tested. One of these models should serve as prototype and the 
other one as scale model. The experimental data of the model should be used to predict 
the behavior of the other one. This increases the confidence factor. 


CONCLUSIONS AND RECOMMENDATIONS 

While several experimental studies can be done on prototype or large scale models, 
size constraints will limit experimentation on complete structures in most cases to small 
scale models. Small scale models can be used as a complement for analytical and 
computational investigations in solving the design problems of complex structures. The 
investigation in this study indicates that the use of small scale models can predict the 
behavior of the prototype very well. 

Partial similarity based on direct use of governing equations is more convenient than 
dimensional analysis, because additional relationships are not needed and the derived 
similarity conditions are based on satisfaction of the field equations of the system. 

In the present study, similitude theory was used by employing systems for which the 
experimental results are not known. In this case, one system was considered to be the 
prototype and another its scale model. Then through the use of the proper scaling laws the 
theoretical data of the model were used to predict the behavior of the prototype. Success 
was measured by comparing the predicted behavior to the analytical results. 

Some recommendations for future research: 

• develop the method for designing and employing scale models for more complex 
systems, i.e. stiffened and/or laminated curved configurations; 

• one of the major problems associated with inelastic analysis of the small scale 
model is the effect of size. A need exists to evaluate the size effect in material 
behavior expecially for geometries with higher scale factors. This is also true in 
dealing with establishing the strength of laminated structures, since strength is 
affected by the accumulation of damage. 
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Cross-Ply Laminated Plates 
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Abstract. The increasing use of laminated composite components for a wide variety of 
applications in aerospace, mechanical and other branches of engineering requires extensive ex- 
perimental evaluation of any new design. Thus, it is extremely useful if a full-scale structure can 
be replaced by a similar scaled-down model which is much easier to work with. The objective of 
this study is to investigative problems associated with design of scaled models. Similitude theory 
is employed to develop the necessary similarity conditions. Both complete and partial similarity 
are discussed. The procedure consists of systematically observing the effect of each parameter 
and corresponding scaling laws. Then acceptable intervals and limitations for these parameters 
and scaling laws are discussed. In each case, a set of valid scaling factors and corresponding re- 
sponse scaling laws that accurately predict the response of prototypes from experimental models 
is introduced. Particular emphasis is placed on the cases of buckling of rectangular cross-ply 
laminated plates under uniaxial compressive and shear loads. This analytical study indicates 
that distorted models with a different number of layers, material properties, and geometries than 
those of the prototype can predict the behavior of the prototype with good accuracy. 



Nomenclature 


a 

plate length 

Q ij > Qij 

lamina stiffness elements 

Au 

laminate extensional stiffnesses 

R 

aspect ratio 

b 

plate width 

u,v,w 

reference surface displacements 

B tJ 

laminate coupling stiffnesses 

t 

ply thickness 

Da 

laminate flexural stiffnesses 

A, 

scale factors 

Ea 

Young’s moduli of elasticity 

Vi : 

Poisson’s ratios 

h 

total laminate thickness 

m 

model 

K t , K s 

non-dimensional critical loads 

P 

prototype 

M x ,M y 

moment resultants 

pr. 

predicted 

N XXl Nyy 

Nry 

inplane normal loads 
inplane shear load 

th. 

theoretical 


INTRODUCTION 


The last step in the design process, before going to production, is the verification of 

t Professor and Head, ’ Graduate Research Assistant, Department of Aerospace Engineering and 
Engineering Mechanics, University of Cincinnati, Cincinnati, OH 45221. 


265 


PREOiOfNG PAGE BLANK NOT FILMED 



266 


DESIGN 


the design. This step necessitates the production of large components and full scale pro- 
totypes in order to test component and system analytical predictions and verify strength 
and performance requirements under the worst loading conditions that the system is ex- 
pected to encounter in service. A scaled-down (by a large factor) model, scale model , 
which closely represents the structural behavior of the full-scale system, prototype , can 
prove to be an extremely beneficial tool. This possible development must be based on 
the existence of certain structural parameters that control the behavior of the structural 
system when acted upon by static and/or dynamic loads. If such structural parameters 
exist, a scaled-down replica can be built, which will duplicate the response of the full-scale 
system. The two systems are then said to be structurally similar. The term, then, that 
best describes this similarity is Structural Similitude. 

Due to special characteristics of advanced reinforced composite materials, they have 
been used extensively in weight efficient aerospace structures. Since reinforced composite 
components require efficiency and wisdom in design, sophistication and accuracy in anal- 
ysis, and numerous and careful experimental evaluations, there is a growing interest in 
small scale model testing[4]. 

By applying similitude theory, we try to find a set of conditions between two similar 
structural systems (scaling laws). Later, these conditions can be used to design a model, 
the experimental data of which can be projected in order to predict the behavior of the 
prototype. 

The objectives of the investigation described herein are: 

• create necessary similarity conditions in order to design an accurate distorted model 

— distortion in stacking sequence and number of plies (N) 

* ply - level scaling 

* sublaminate - level scaling 

— distortion in material properties £,j, v HiP 

• evaluate the derived similarity conditions analytically. 

Similarity conditions provide the relationship between model and its prototype, and 
can be used to extrapolate the experimental data of a small and less expensive model in 
order to predict the behavior of the prototype. This study presents the applicability of 
small scale models, especially distorted models, in analyzing the elastic behavior of cross- 
ply laminated plates. Furthermore, it is assumed that the laminates are free of damage 
(delaminations, matrix cracking, fiber breaks, etc.). 

In this study, we consider only the procedure that is based on the direct use of the 
“governing equations”. This method is more convenient than dimensional analysis, since 
the resulting similarity conditions are more specific and the relationships among variables 
are forced by the governing equations of the system. 

Often complete similarity is difficult to achieve or even undesirable. When at least 
one of the similarity conditions can not be satisfied, partial similarity is achieved. In this 
case, the model which has some relaxation in similarity conditions is called a distorted 
model. Distorted models are more practical, since relaxation of each similarity condition 
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eliminates some restrictions on the model design. These relaxations in the relationship 
between two systems cause model behavior to be different from that of the pro o ype. 
Since each variable has different influence on the response of the system, the resu mg 
similarity conditions have different influence. By understanding the effect of vana es 
and similarity conditions over desired intervals, the similarity conditions which have the 
least influence can be neglected without introducing significant error[3j. 

Buckling of Symmetric Laminated Cross-Ply Rectangular Plates 

Consider symmetric cross-ply laminated plates {B i: = 0 ,D 16 = D ™ir 
0) The plates are subjected to inplane normal and shear loads ( iV xx , iv yv , iv xy j. e 

governing differential equations for buckling and vibration of symmetric cross-p y rectan- 
gular plates are as follow [1] 

^ll u !Lr + A 12 V° yx + A66(u° yy + u ?ry) = 0 ^ 

+ v% t ) + Anu 0 ^ + A 22 v° yy = 0 (2) 

D\\W° XXXX + 2^i2^°irvv ^ 22W °yyyy ~ ^ xx w * x ~~ ^yy w wv ~ ^ x v w ,xy ~~ ^ 

where D\ 2 = P 12 -b 2Z^66 

For simply supported plate, the approximate boundary conditions are 


at x = 0, a 


= 0 , M x = —Dn w % x = 0 


at y = 0, b w — 0 , , M y — D 22 w yy 0 

J Uniaxial Load : Consider the plates to be subjected to inplane uniaxial 

compression load in x direction ( N xx ). The buckling differential equation is: 

+ 2fi 12 to°„ vv + D !2 w°„ wy - N„w% = 0 (4) 


The solution 


. mu . ,mry 

w = A mn szn( )sm(— 7 — ) 

a 0 


satisfies the boundary conditions. Then, 


N„ = ADu(t? + 2^4? + D*(t) 4 0’1 


By applying similitude theory to Eq.( 5) 


. _ . At _ A> _ A& 

A 2 - Dv X 4 X ? m 


which yields to following scaling laws 




■^£22 K 


( 7 ) 
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Ar-„ = Xi-A^i 


AeV* 


A K.. = 

K 


A^A 3 '" 1 R 


(S) 

( 9 ) 


where K TX = 


N*:b 2 
E 22 h* 


Boundary conditions do not give any similarity condition. These conditions, Eqs.( 7) — 
( 9), involve response^*-,,, A m , A n ) and structural geometric parameters (A Dij , X h , X R ). 

H '• Shear Buckliri£J We now consider a simply supported plate which is sub- 

jected to the an in-plane shear stress ( N xx = N yy = < = 0). 


Dn w °x xxx + 2£h2U>° zyv + D 22 w° yyyy — 2N xy w° xy = 0 


(10) 


The solution of the form 


w 


o v-' «. . . .mil, . ,nxy x 

= n ^mn5in( )szn(— — — ) 

m=ln=l d ° 


satisfies all B.C.’s, but does not satisfy the buckling equation, Eq.( 10). Use of the 
Galerkin procedure, yields 


r _Dn_m< ( 0 J> 12 mV , D 22 D 4 , , 32mn ^ ^ „ 

( E 22 h 3 it! 3 + 2 E 22 h 3 R + E 22 h 3 Rn ' Amn ~ Kl ^ £ A PiQ™ Pq 

p=l 7=1 

for m, n = 1, 2, • • • , co ; subject to the constraints m ± p = odd and n ± q = odd. 
v - ^v &2 d _ a n M 

3 Enh 3 ’ 6 ’ ( m 2_ p2 )( n 2_ ?2) 

Applying similitude theory to Eq.( 11) 

-A On A« A fl „ Aj^A 3 A^ 


( 11 ) 


ArAI = 




Ae„A 3 A 3 A^A 3 A fl A^A J' w *» 


Am A n A(j 


( 12 ) 


where 


n = ££v? 

p=l ? =l 


mnp<j 


Parenthesis: 


m,n,p, and are integers which depend on the number of terms needed to approxi- 
mate well the buckling mode shape (Symmetric/ Anti-symmetric). By assuming the same 
aspect ratio for the model and its prototype (A^ =1) and similar construction, model 
and prototype both can be well approximated by the same number of terms in the series 
with the same contribution of terms to the buckling mode and thus 


Am — A n — \ — \ Q — 1 


An 1 , A^ — 1 
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Eq.( 12) yields to the following scaling laws 


A*. = 
A k, = 
A k. = 


K Du 


A£ m A^Ar 

Agn 

A e 2 j A^ An 

A Dr: \ 


(13) 

(14) 

(15) 


These 

eters. 


conditions, Eqs.( 13)-( 15), also involve response and structural geometric param- 


Comvlete Si milarity : 

The necessary condition for complete similarity between the model an i s pr yp 
is that all scaling laws be satisfied. This requirement yields 

i _ \ _ i (16) 

ADu — Ap,2 — *D?2 

Bv inspection it can be seen that these conditions, Eqs.( 16), are independent of ply 
thlkness and they only depend on materia! properties and TO /g P * “ 

with different ply thickness but the same stacking sequences i.e.(0/90). and (U./9U.M 
a isfv Eos ( 16) This is called ply-level scaling and it is the easiest way to achieved 
complete similarity. Table 1 presents the ply-level scaling for cross-ply plates under inplane 

^lo'fatwe have shown that in the special case of ply-level scaling similarity can be 
achieved. However, there are some constraints, in designing the mo e ese cons r un 
involve the geometry of the model, the model material, the number of plies and the 
stacking sequence of laminates. Since this still appears to be restrictive, we allow the use 
of distortion in the design of the model. 

Table 1 Comparison of shear hackling loads of Kevlar/Epoxy plates with ply - level 

scaling (complete similarity) . 


model 

_ NzyV 

%D 

isc. 

Configuration 

model 

prototype 

predicted 

th.{p)kpr.{p) 

th.(p)&th.(m) 

(0o/902)5 

32.74 

32.74 

32.74 

0.0 1 

0.0 

(0io/90io)j 

32.74 

32.74 

32.74 

0.0 

0.0 

(02O/902O) J 

32.74 

32.74 

32.74 

0.0 

0.0 


| theory — predicted | 

%Disc.{th.kpr.) = 100 x ^ 

Partial Similarity : When at least one of the design scaling laws cannot be 

satiS^Tp^l similarity is achieved. By understanding the > effect of “ d 

accuracy of the scaling laws over desired intervals, the scaling aws w ic yie 
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accurate prediction for the prototype can be chosen as corresponding to the “right” type 
of distortion. 

The choice of the “right” type of distortion is investigated as follows. In each case, all 
of the model parameters except one, are chosen to be identical to its prototype. Then, 
the effect of this relaxation for a wide range of this parameter is investigated. 

Nuniber of Plies. / There are three ways to scale down the number of plies in 
a model, a) ply-level scaling((0 n /90 n ) J ) b) sublaminate level scaling[2]((0/90) nj ) and 
c) general reduction of plies. The ply-level scaling leads to complete similarity(as al- 
ready discussed). But the two other methods yield partial similarity. Figures 1 and 
2 present the compressive and shear buckling loads for models with different number of 
plies than those of the prototype and the predicted loads by using scaling laws depicted by 
Eqs.( 8 )and( 14). Both sublaminate level and general scaling are presented. It is shown 
that all models(except for (0/90/0)) can predict the prototype behavior with excellent 
accuracy. 

Table 2 Accuracy of models with sublaminate level scaling (;V p ^ N m ) using 
non-dimensional load I\ } , \e, } = X u>j = 1; (0/90)2Oj- 


model 

K - ML 

' Ellh 3 

%Disc . 

Configuration 

model 

prototype 

predicted 

th.(p)kpr.(p) 

th.(p)kth.(m) 

(0/90)2, 

32.740 

34.040 

32.740 

3.82 

3.82 

(0/90)5, 

34.009 

34.040 

34.009 

0.09 

0.09 

(0/90)i O j 

34.034 

34.040 

34.034 

0.02 

0.02 


Table 3 Accuracy of models with sublaminate level scaling (N p N m ), using 
non-dimensional load I\,, Xg tj = \ Vij = 1 ; (0/90)2o 4 - 


model 

jt- - N ^‘ 
1 D72 

%Disc. 

Configuration 

model 

prototype 

predicted 

th.(p)kpr.{p) 

th.(p)kth.(m) 

(0/90)2, 

150.20 

56.763 

54.594 

3.82 

164.61 

(0/90)5, 

63.045 

56.763 

56.709 

0.10 

11.07 

(0/90)io, 

58.718 

56.763 

56.752 

0.02 

3.45 


Material : Now we consider distortion in model material. For this purpose two 
different groups are considered: Isotropic materials (which include metals and plastics), 
and fiber reinforced composites. In all of these cases the prototype is an orthotropic 
laminated plate. 

For the composite model, model and prototype have the same stacking sequence, number 
of plies ((0/90/0. -Ou) and aspect ratio. Figures 3 and 4 present theoretical and predicted 
buckling loads of prototype and theoretical ones of the models for some typical composite 
materials. For the Kevlar/Epoxy prototype a Boron/Epoxy, Boron/ Polymide, and most 
of Graphite/Epoxy’s can be used as the model material or vice versa. But Glass/Epoxy 
is not a good choice for predicting a Kevlar/Epoxy prototype or vice versa. 
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Since plastics are used extensively in experimental studies of the behavior of the struc- 
tures, the possibility of a plastic model or in general a model with isotropic materials is 
considered. For isotropic materials, the assumption of = 1 yields a model which cannot 
predict accurately the behavior of the prototype. By choosing R m as a design parameter 
we are able to find isotropic models which yields excellent accuracy. Scaling laws depicted 
by Eq.( 8) and Eq.( 14), yield acceptable aspect ratios for the models. Tables 4 and 5 
present theoretical and predicted buckling loads for prototypes when the corresponding 
models are isotropic material. 

Table 4 Relaxation in material properties by using isotropic model. 

Prototype Is Kevlar/Epoxy (0/90)2o» ■ 


model 


%Disc. 

material 

Rm 

model 

prototype 

predicted 

th.{p)kpr.{p) 

th.(p)kth.(m) 

Aluminum 

0.705 

4.174 

14.16 

14.16 

0.0 

70.51 

Brass 

0.705 

4.143 

14.16 


0.0 

50.80 

Copper 

0.705 

4.149 


MI1MI 

0.0 

70.69 

Steel 

0.705 

4.045 

14.16 

mam 

0.0 

71.36 

PVC 

0.705 


14.16 

14.16 

0.0 

69.1 

Polyethylene 

0.705 

4.760 

14.16 

14.16 

0.0 

66.38 


Table 5 Relaxation in material properties by using isotropic model. Prototype is 

Kevlar/Epoxy ( 0 / 90 ) 2 Oj- 


model 



%D isc. 

material 

Rm 

model 

prototype 


th.{p)kpr.{p) 

th.{p)k.th.{m) 

Aluminum 

0.627 

20.22 

34.04 

34.04 

0.0 

40.61 

Brass 

0.627 1 

20.07 

34.04 

McgligM 

0.0 

41.02 

Copper 

0.627 1 

20.10 

34.04 

34.04 

0.0 

40.97 

Steel 

0.627 

19.63 

34.04 

34.04 

0.0 

42.32 

PVC 

0.627 


34.04 

34.04 

0.0 

37.77 

Polyethylene 

0.627 | 

23.05 

34.04 

34.04 

0.0 

32.28 


DISCUSSION 

An analytical investigation has been conducted in order to establish the applicability of 
similitude theory to cross-ply laminated plates. The results presented herein indicate 
that for elastic response of a cross-ply rectangular plate, based on structural similitude, 
a set of scaling laws can be found to develop design rules for small scale models. By 
establishing similarity conditions, the model parameters are specified. The accuracy of 
predicted prototype behavior by various models is investigated. The verification is based 
on the exact analytical solution of the mode! and its prototype. Theoretical compressive 
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DESIGN 


and shear buckling loads of scale models are projected to prediet corresponding 
behavior. Even for models with different number of plies and stacking sequences ( 
models) the predicted data are well matched with full scale prototype data, or bo h load 
conditions the enforcement of the same response scaling aw q.( ) an f 

models which yield accurate predictions. In the present study we assume that ex pt or 
isotropic models, prototype and models have the same aspect ratio(A R - l). 1 his n 
necessary condition and in general (A, / 1). The accuracy of prediction ,s very sensitive 

to the scale factor of the aspect ratios. 

CONCLUSIONS AND RECOMMENDATIONS 

Small scale models can be used as a complement for analytical and computational investi- 
gations in solving the design problems of complex structures. This study mtoto i that a 
distorted model with a fewer number of layers can predict buckling load of the prototype 

W ‘ l pfrt°al stoilarity based on direct use of governing equations is very convenient. There 
is tremendous freedom in design scale models because the number of simi arity con i 10 
is much smaller than the number of design variables. 

Some recommendations for future research. 

• Develop the method for designing and employing scale models for more complex 
systems, i.e stiffened and/or laminated curved configurations. 

• Experimental verification of the accuracy of the purposed scaled model. 

• Implements of the structural similitude to inelastic and failure analysis of com- 
posite structures. 
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Figure 1: Predicted and theoretical compressive buckling load of the prototype (0/90) 
when (0/90/0. ..) n is used as model. 
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Figure 2: Predicted and theoretical shear buckling load of the prototype (0/90) 
(0/90/0. ..) n is used as model. 








Figure 3: Predicted and theoretical compressive buckling load of the Kevlar/Epoxy 
prototype when model have different material properties. 



Figure 4: Predicted and theoretical shear buckling load of the Kevlar/Epoxy prototype when 
model have different material properties. 


274 



DESIGN OF SCALED DOWN MODELS FOR 
STABILITY AND VIBRATION STUDIES 


J. Rezaeepazhand * and G. J. Simitses 

Department of Aerospace Engineering and Engineering Mechanics 
University of Cincinnati, Cincinnati, OH 45221. 

ABSTRACT 

Use of reinforced composites in light-weight aerospace structures has increased stadily over 
the years. The outstanding mechanical and physical properties of advanced composites 
provide the engineer with potential to optimize properties specific to application. Since 
reinforced composite components require efficiency and wisdom in design, sophistication 
and accuracy in analysis, and numerous and careful experimental evaluations, there is a 
growing interest in small scale model testing 1 . 

A scaled-down (by a large factor) model, scale model , which closely represents the struc- 
tural behavior of the full-scale system, prototype , can prove to be an extremely beneficial tool. 
This possible development must be based on the existence of certain structural parameters 
that control the behavior of the structural system when acted upon by static and/or dynamic 
loads. If such structural parameters exist, a scaled-down replica can be built, which will du- 
plicate the response of the full-scale system. The two systems are then said to be structurally 
similar. The term, then, that best describes this similarity is Structural Similitude. 

Similitude theory is employed to develop the necessary similarity conditions (scaling laws). 
Scaling laws provide relationship between a full-scale structure and its scale models, and can 
be used to extrapolate the experimental data of a small, inexpensive, and testable model into 
design information for a large prototype. The difficulty of making completely similar scale 
models often leads to accept certain type of distortion from exact duplication of the prototype 
(partial similarity). Both complete and partial similarity are discussed. The procedure 
consists of systematically observing the effect of each parameter and corresponding scaling 
laws. Then acceptable intervals and limitations for these parameters and scaling laws are 
discussed. In each case, a set of valid scaling factors and corresponding response scaling laws 
that accurately predict the response of prototypes from experimental models is introduced. 

* Graduate Research Assistant. * Professor and Head, Associate Fellow of AIAA. 
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Particular emphasis is placed on the cases of free vibration and buckling of rectangular 
angle-ply laminated plates under uniaxial compressive and shear loads. This analytical 
study indicates that distorted models with a different number of layers, material properties, 
and geometries than those of the prototype can predict the behavior of the prototype with 
good accuracy. 

The objectives of the investigation described herein are: 

• create necessary similarity conditions in order to design a distorted model that accu- 
rately predicts prototype behavior. 

— distortion in stacking sequence and number of plies (N) 
ply - level scaling and sublaminate - level scaling 

— distortion in material properties Eij, Vij,p 

— distortion in fiber orientation angle 9 


• evaluate the derived similarity conditions analytically. 


In all of our work in this area we will restrict ourselves to linearly elastic material behavior. 
Furthermore, it is assumed that the laminates are free of damage (delaminations, matrix 
cracking, fiber breaks, etc.). 

In this study, we consider only the procedure that is based on the direct use of the 
“governing equations”. This method is more convenient than dimensional analysis, since the 
resulting similarity conditions are more specific and the relationships among variables are 
forced by the governing equations of the system. 

Consider symmetric angle-ply laminated plates (B,j = 0). The plates are subjected 
to inplane normal and shear loads ( N„ , N yy , N x y ). The buckling loads and vibration 
frequencies of symmetric angle- ply rectangular plates axe described only by one differential 
equation 2 : 


D xl w Q „ SI +AD l& w 


o 

,xxxy 


+ 2f} 12 ty° xxyv +4l)26^° C yy V + -D22« j0 11 


,yyyv 


-NnW^-NyyW^-N^W 


'rXX 


’ xy ^ r ry 


0 

= P W ,tt 

(i) 


where Du = Du + 2-D&6 


For a simply supported plate, the boundary conditions are 

at x = 0, a w = 0 , M x = —Dn^xx = 0 


2 



at y = 0, b w — 0 , , M y — — Z^ju; 0 ^ — 0 

I. Shear Buckling : We now consider a simply supported plate which is subjected 
to the an in-plane shear stress ( N ; ** = N yy = w° tt = 0). 

DuW° xxxx + 4 D\sUp xxxy + 2Di2W° xxyy + 4£>26 u > r xyj,y + DK W ,yyvy ~ ~ ® ^ 

The solution of the form w° = £“ =1 A mn sin(^)sin(^) satisfies all B.C.’s, but does 
not satisfy the buckling equation, Eq.( 2). Use of the Galerkin procedure, yields 


m’ - m 2 n 2 _ n’. a 

(U*ii— T + + D 22 —)A mn - 


^ £ £( (m; t — o,6 + 

ir 2 ab ±!,tl a 2 


-^2 ~-^26 + 


(3) 


for m, n = 1, 2, • • • , co ; subject to the constraints m ± p = odd and n ± q - odd. 




where R=-^ , Qm npq - (m , _ ^,,2 _ g2) 

m, n, p, and q are integers which depend on the number of terms needed to approximate 

well the buckling mode shape(Symmetric/Anti-symmetric). By assuming the same aspect 
ratio for the model and its prototype (A* = 1) and similar construction, model and prototype 
both r*.n be well approximated by the same number of terms in the series with the same 
contribution of terms to the buckling mode and thus 


Am — A n A j> A q 1 


Aft = 1 J ^A mn ~ 1 


Applying similitude theory to Eq.( 3) yields the following scaling laws 

Ad u 


A/r. = 
A k. = 
A/c. = 
A k. = 


Api, 

AsmA^Ah 
A Dm 


Adj 2 A^ 

Ad i» 


A/l 


A/C. = 


A/^A 2 r 

Ad,» 


AejjA^ 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


where K, = 




£ 22 /i 3 
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These conditions, Eqs.( 4)-( 8), contain both response and structural geometric param- 
eters. 

Complete Similarity : The necessary conditions for complete similarity between 

the model and its prototype is that all seeding laws be satisfied. This requirement yields 

By inspection it can be seen that these conditions, Eqs.( 9), are independent of ply thicknesses 
and. they only depend on material properties and number of plies. So, two plates with 
different ply thickness but the same stacking sequences (i.e.(0/ + 45/90/ - 45), and (0„/ + 
45«/90 n / - 45„),) satisfy Eqs.( 9). This is called ply-level scaling and it is the easiest way to 
achieved complete similarity. Table 1 presents the ply-level scaling for angle-ply plates under 
inplane shear loads. 

So far, we have shown that in the special case of plv-level scaling similarity can be 
achieved. However, there are some constraints, in designing the model. These constraints 
involve the geometry of the model, the model material, the number of plies and the stacking 
sequence of laminates. Since this still appears to be restrictive, we allow the use of distortion 
in the design of the model. 


Table 1 Comparison of shear buckling loads of Graphite/Epoxy plates with ply - level 

scaling (complete similarity). 


model 

v Si 

A> - Erih 3 . 

%D tsc. 

Configuration 

model 

prototype 

predicted 

th.(p)kpr.(p) 


(0/ + 45/90/ - 45), 

30.721 

30.721 

30.721 

0.0 

0.0 

(0 3 / -1- 45 3 /90 3 / - 45 3 ), 

30.721 

30.721 1 

30.721 

0.0 

0.0 

(0io/ + 45io/90io/ — 45io)j 

30.721 

30.721 

30.721 

0.0 

0.0 


where %Di<c.{th.!v ,.) = 100 x 


Partial Similarity : Often complete similarity is difficult to achieve or even unde- 

sirable. When at least one of the similarity conditions can not be satisfied, partial similarity 
is achieved. In this case, the model which has some relaxation in similarity conditions is 
called a distorted model. These relaxations in the relationship between two systems cause 
model behavior to be different from that of the prototype. Since each variable has different 
influence on the response of the system, the resulting similarity conditions have different 


4 













influence. By understanding the effect of variables and similarity conditions over desired 
intervals, the similarity conditions which have the least influence can be neglected without 
introducing significant error 3 . 

The choice of the “right" type of distortion is investigated as follows. In each case, all 
of the model parameters except one, are chosen to be identical to its prototype. Then, the 
effect of this relaxation for a wide range of this parameter is investigated. 

Number of Plies : There are three ways to scale down the number of plies in a model, 
[a] ply-level scaling (0 n /+45 n /90 n /-45 n ),) [bj sublaminate level scaling 4 ( (0/+45/90/-45) n ,) and 
[c] general reduction of plies. The ply-level scaling leads to complete similarity(as already 
dismissed). But the two other methods yield partial similarity. Figure 1 presents the shear 
buckling loads for models with different number of plies than those of the prototype and the 
predicted loads by using scaling laws depicted by Eq.{ 6). 


Table 2 Accuracy of models with sublaminate level scaling [Np A Nm) using 
non-dimensional load K s , A e,, — = 1 / ( 0 / + 45/90/ — 45)ioj- 


model 

u 

II 

Cr 

CO ►- 

%L 

)i,sc. 

Configuration 

model 

prototype 

predicted 

th.{p)kpr.(p) 

th.{p)kth.(m) 

fo/ + 45/90/ - 45)2, 

37.813 

44.218 

44.714 

1.12 

14.48 

(0/ +45/90/ - 45) 3 , 

40.391 

44.218 

44.187 

0.25 

5.63 


Material : Now we consider distortion in model material. For this purpose two 
different groups are considered: Isotropic materials (which include metals and plastics), 
and fiber reinforced composites. In all of these cases the prototype is an angle-ply laminated 

plate. 

For the composite model, model and prototype have the same stacking sequence, number 
of plies ((0/ + 45/90/ - 45) i 0 .) and aspect ratio. Figure 2 presents theoretical and predicted 
buckling loads of prototype and theoretical ones of the models for some typical composite 
materials. For the Kevlar/Epoxy prototype almost all considered materials can be used as 
the model material or vice versa. 

Since plastics are used extensively in experimental studies of the behavior of the struc 
tures, the possibility of a plastic model or in general a model with isotropic material is 
considered. For isotropic materials, the assumption of A* = 1 yields a model which cannot 
predict accurately the behavior of the prototype. By choosing Hm as a design parameter 
we are able to find isotropic models which yields excellent accuracy. Scaling law depicted 
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by Eq.( 6), yields acceptable aspect ratios for the models. Tables 3 present theoretical and 
predicted budding loads for prototypes when the corresponding models are made of isotropic 

material. 

Table 3 Relaxation in material properties by using isotropic model, 


prototype is Kevlar/Epoxy (0/ + 45/90/ 45)io* • 


model 

N-,b~ 

A * " E~f i 3 

%L 

isc . 

material 

Rvn I 

model 

prototype j predicted 

th.(p)kpr.(p) 

th.(p)kthjm) 

Aluminum 

0.868]' 

12.08 

43.34 ] 43.34 

o.o ~ 1 

72.12 

Copper 

0.868 

12.01 

43.34 | 43.34 1 

0.0 

72.29 

PVC 

0.868 

12.66 

43.34 | 43.34 1 

0.0 

70.79 


This study presents the applicability of small scale models, especially distorted models, 
in analyzing the elastic behavior of angle-ply laminated plates. Distorted models axe more 
practical, since relaxation of each similarity condition eliminates some restrictions on the 
model design. The results presented herein indicate that, for elastic response of an angle-ply 
rectangular plate, based on structural similitude, a set of scaling laws can be found to develop 
design'rules for small scale models. Results for buckling characteristics under uniaxial and 
vibration characteristic will be presented in the full paper. 
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N m (Number of Plies) 


Figure 1: Predicted and theoretical shear buckling load of the prototype (0/ + 45/90/ 45)io, 

when (0/ +45/90/ — 45) nj is used as model. 



Figure 2: Predicted and theoretical shear buckling load of the Kevlar/Epoxy prototype when 
model have different material properties. 
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